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Abstract 



We analyze accumulation of ultracold neutrons (UCN) in a supcrfluid-hclium converter vessel 
surrounded by a magnetic multipole reflector. Employing formulas valid for trapped UCN in 
mechanical equilibrium we solved the rate equation and obtained results for the saturation 
density of low field seeking UCN for any field strength of the multipolar field. The addition 
of magnetic storage can increase the density and energy of the UCN produced and serves to 
mitigate the effects of wall losses on the source performance. This work was performed in 
preparation of the UCN source project SuperSUN at the ILL. 
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1 Introduction 

■^j- ■ Mirror reflection of neutrons is an effect of the neutron optical potential which is mainly due 

to coherent s-wave scattering of neutrons by atomic nuclei in condensed matter pQ. Neutrons 
with energy sufficiently low to become totally reflected under any angle of incidence are called 
ultracold neutrons (UCN). This peculiar property enables experimentalists to store them in 
"neutron bottles" made of suitable materials with small cross sections for neutron absorption 
and providing well depths up to about 300 neV [21 [3]. Storage time constants of many hundreds 
of seconds and the possibility to employ also magnetic fields and gravity for trapping and 
manipulation have made UCN a versatile tool to investigate various phenomena of fundamental 
physics complementary to experiments at high-energy particle accelerators [HO [6]. 



Among recent experiments with UCN feature a first demonstration of gravity resonance 
spectroscopy with the goal to search for deviations from Newton's gravity law at the micrometer 
length scale [7], searches for "mirror dark matter" [8, 9j, a test of Lorentz invariance [10], searches 
for axion-like particles [11[ \12\ [T3] , a demonstration of the effect of accelerated matter on the 
neutron wave [14] and of the stability of the Berry geometrical phase for spin ^ particles under 
the influence of noise fluctuations [15]. Earlier work with UCN on the geometrical phase was 
published in |17[ 118], while its first demonstration with cold neutrons can be found in [16] . 

Long standing are efforts to improve the accuracy of the weak axial- vector and vector coupling 
constants of the nucleon derived from precise values of the neutron lifetime [23 [2D [221 [23] 
and the beta-asymmetry [241 [25| [26] [27] . Among other applications these values are crucial input 
for calculations of weak reaction rates in big-bang nucleo-synthesis and stellar fusion [28, 29], and 
of the efficiency of neutrino detectors [30J . Also long standing is the search for a non- vanishing 
neutron electric dipole moment (EDM) , which would violate the symmetries of time reversal (T) 



1 



and thus via the CPT theorem also the combined symmetries of charge conjugation and parity 
(CP). This search was proposed already in 1950 by Purcell and Ramsey [31] and has become 
a prominent route to investigate new mechanisms of CP violation beyond the standard model 
complex phase of the weak quark mixing CKM matrix, and the matter-antimatter asymmetry in 
the universe [32] . At the present best level of sensitivity, still no EDM was observed [33] . Several 
projects are in preparation or underway with the goal to gain at least an order of magnitude in 
sensitivity [3H [351 ESJ EH [381 EHl SQl SI] . 

Most studies with UCN are counting statistics limited and will strongly profit from new 
UCN sources which are currently being developed in various laboratories around the world 
g2l S3l SH SSI S71 SHI S9]. They are all based on the " superthermal" UCN production 
scheme proposed in 1975 by one of the authors together with Mike Pendlebury [50] . using either 
superfluid 4 He or solid deuterium as a medium for neutron conversion. Here we focus on UCN 
production in a converter of superfluid 4 He installed at the end of a neutron guide, wherein 
neutrons with energy 1 meV, respectively, wavelength 0.89 nm may loose nearly their entire 
energy in single scattering events. At low temperature only few excitations are present in the 
helium that are able to scatter UCN back to higher energies. With the vanishing absorption 
cross section of 4 He it becomes possible to accumulate UCN within a converter with reflective 
boundaries before releasing them to an experiment at room temperature. While an earlier 
attempt to bring this technique to life was hampered by extraction losses (nonetheless producing 
record UCN densities for its time) [51], a more efficient method was developed recently by one 
of the authors together with his co-workers [491 l52j [53] . UCN are extracted through a cold UCN 
valve and a short vertical UCN guide section, superseding lossy separation window, screens and 
gaps for thermal insulation between the converter and the UCN guide of the earlier scheme. 
Work is in progress to bring the technique to maturity for a new user facility at the ILL, and 
in particular to perform a neutron lifetime experiment using magnetic trapping [541 [55| 156] . 
Also other groups have recognized the potential advantages of a superfluid helium converter 
feeding UCN to experiments at room temperature [431 147j. and in some experiments this type 
of converter is employed in situ [311 ESI [57] . 

The efficiency of a UCN accumulator at an external neutron beam relies on loss rates being 
as low as possible. Most critical are those losses which occur when UCN hit the walls of the 
converter vessel. They are proportional to the frequency of collisions and thus depend on the 
size and shape of the converter vessel. From transmission measurements through superfluid 4 He 
at 1.25 K a mean free path of 17 m was derived for the 0.89 nm neutrons most effective for 
UCN production [58J. Hence, the vessel can be made several meters long without significant 
reduction in UCN density. On the other hand, the lateral dimensions of the converter vessel 
should match the size of the available neutron beam and guide it to avoid dilution of the incident 
flux. The mean free path of UCN in a vessel with such geometry is therefore only in the order 
of 5 — 10 cm, leading to high typical frequencies of UCN wall collisions of 50 — 100 per second. 
It thus becomes challenging to obtain long UCN storage time constants which however are 
a prerequisite for accumulation of a high saturated UCN density. Values measured for narrow 
vessels are normally well below 200 s. For instance, in a recent experiment on UCN production, a 
rather short storage time constant of 67 s was obtained for a vessel held at 0.7 K, which consisted 
of a 1 m long 7x7 cm 2 tube of BeO with Be windows on each end and included a short pipe 
from stainless steel. That, nonetheless, a record UCN density was obtained demonstrates the 
potential of the method [49J. To our knowledge the Cryo-EDM collaboration achieved with 
t = 160 s the so far highest value for a helium converter enclosed within matter boundaries, 
using a 3 m long tubular vessel with diameter 63 mm, made of Be coated copper and closed off 
by Be windows [34] . 

Magnetic trapping of UCN offers a viable way for a drastic improvement of the storage prop- 
erties of the converter vessel, ultimately limited only by the neutron lifetime ~ 880 seconds. 
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It relies on the potential energy ±pB of the neutron magnetic moment fi of 60 neV/T in a 
magnetic field B. Suitable configurations of magnetic field gradients keep low field seeking UCN 
away from walls where otherwise the collisional losses occur. A group at NIST has demonstrated 
storage time constants consistent with the neutron lifetime within a helium converter equipped 
with a superconducting magnetic quadrupole UCN reflector [571 [59]. The apparatus was de- 
signed to perform neutron lifetime measurements, for which a complete suppression of UCN wall 
contacts is mandatory. On the other hand, for the sake of enhancing the output of a UCN source, 
combined magnetic and material trapping turns out to be particularly beneficial. In addition, 
the phase space for UCN accumulation can be much increased using a higher multipole order. 

In this paper we provide an analytic treatment of the rate equation for UCN production and 
storage in a superfluid-helium converter confined by material walls and surrounded by a magnetic 
mirror. We show that, combining a converter vessel possessing good (but not exceptional) storage 
properties with a magnetic mirror of high multipole order, one may achieve a saturated UCN 
density close to the theoretical limit defined by an ideal experimental bottle, i.e. a square well 
potential without imaginary part. That the magnet needs to generate only part of the trapping 
potential is of great practical value for constructing a device using standard superconducting 
wire technology. 



2 Rate equation and system definition 

The temporal evolution of the spectral UCN density in a closed converter irradiated with the 
cold beam is governed by a simple rate equation. UCN production is characterized by a spectral 
rate density p that depends on the spectral flux of the incident beam, and a loss term is due to 
finite lifetime r of UCN in the converter, 

dn(e ,i) , . n(e ,t) 

Here we label stored neutrons by their total energy eo, defined as their kinetic energy at the 
point of lowest potential energy in the trap. The quantities n(eo,t) and p(eo) denote the real 
space density, respectively, the spatial UCN production rate density, each per energy interval 
of a group of UCN with total energy in the range (eo, eo + deo). The saturated spectral UCN 
density obtains when UCN losses balance UCN production for t ^> r after having switched on 
the beam. It is given by 

noo (eo) = p (e ) t (e ) ■ (2) 

If one wants to fill a trap with UCN up to a cutoff energy set by the trapping potential Vt rap , 
what matters is the saturated total UCN density which is obtained from 

/"Vtrap 

noo= (e ) de . (3) 



Jo 

Many applications of UCN sources involve filling external traps with as many UCN as possible, 
followed by a long time for holding or manipulation, during which the density in the source can 
be refreshed. Therefore useful parameter of the converter to be optimized^- 

We consider a system as schematically shown in fig. 1. A cylindrical converter vessel with 
diameter 2R is situated within a multipole magnet and illuminated homogeneously by a cold 
neutron beam passing along the r = axis. UCN accumulation takes place over a length L ^> R 
between a beam window and a UCN valve fully immersed in the helium as in the apparatus 



1 For experiments involving external traps with poor storage properties it will be better to drain UCN frequently 
from a partially charged source. However, also in this case a long UCN storage time constant is an asset as it will 
raise the time-averaged UCN content of the converter. 
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Figure 1: Schematic of the UCN accumulator comprized of a superfluid 4 He converter with 
multipole magnet and system for UCN extraction. On the left a cut view is shown for n = 12; 
filled (open) squares indicate electric current flowing into (out of) the plane. The neutron optical 
potentials are: V for the cylindrical part over the length L, V at the beam window and the 
UCN valve, and > V for the UCN extraction system. 



described in [60 1. Shown is a butterfly valve but also different types may be envisaged, such 
as an iris diaphragm valve. For experiments at room temperature UCN are released into a 
window-less extraction system with a short vertical guide section as described in [530. In the 
section for UCN accumulation the cylindrical wall possesses a neutron optical potential V + iW 
mil, with 

V = —Y,N l b l , W = %Y i N l va l (y), (4) 
mn l 2 i 

where m n is the neutron mass, iVj is the atomic number density of the nuclear species I with 
coherent bound neutron scattering length bi, and 07 (v) is the sum of cross sections for neutron 
capture and scattering at (thermal) neutron velocity v, which is usually proportional to 1/v. 
The beam window and the UCN valve are made of (or coated with) a material with neutron 
optical potential V + iW . 

A radial n-polar magnetic field with modulus 

— — 1 

B (r) = B R (£) 2 (5) 

can be generated as shown with a regular arrangement of an even number of n straight current 
bars on the outer cylinder surface, with opposite current in adjacent bars (in practice one employs 
long racetrack coils). A neutron moving in such a field has a magnetic potential energy of 



r 



f" 1 



V m (r) = ±V mR (-)' , (6) 

where V m R is defined to be positive so that the positive (negative) sign in eq. [6] stands for 
neutrons in the low (high) field seeking spin state. The trapping potential of the converter 
vessel filled with superfluid helium (Vn e = 18.5 neV) is given by 



U tmp = min I V + V mR , V - V Ke . (7) 



2 It is also conceivable to place the UCN valve closer to (or within) the extraction guide. This would however 
considerably increase the surface of wall material exposed to the UCN during accumulation. Here we analyse the 
system as shown in fig. 1. 
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3 UCN losses from the converter 



The inverse of the time constant appearing in the loss term in eq. dj i.e. the rate of UCN loss 
normalized to the number of UCN in the converter, is comprized of several contributions, 



T- 1 



(*o, T) = T~l n (e ) + (e ) + r" 1 (T) + r£ + r^ pol (e ) + rf. (8) 

From left to right, the partial rates are due to UCN loss at collisions with the walls of the 
converter vessel, escape of UCN through apertures caused by manufacturing imperfections of 
the vessel assembly and through slits of an imperfect UCN valve, upscattering by thermal 
excitations in the helium, absorption by 3 He impurities, UCN depolarization at wall collisions, 
and neutron beta decay. Note that the rate 7"^ epo i may become relevant only if trapping is at 
least partly magnetic. In an experimental study a depolarization probability per wall collision 
of 7 x 10~ 6 was measured for a bottle made of Be [61) . This corresponds to a rate smaller than 
r^ 1 even for the fastest neutrons in the narrow trap geometry discussed here. Like the first two 
rates in eq. [8] it will be further suppressed due to the multipole magnet. We therefore neglect 
T depoi- F° r temperatures below 1 K [60| 162). 

*pi-M, (») 



ll|> 1 * 100s 

so that for T < 0.5 K the rate r~p contributes less than 10% of rl 1 . The rate T~ !oa can be 
suppressed below any relevant level by purification of the helium from 3 He using a superleak 
[53J E31 El] or the heat flush technique [65J. As a result, we are left with r~* u and r~ lit as 
dominating contributions, and the rate due to neutron decay sets an ultimate lower limit for a 
perfect converter. 

For the losses due to wall collisions we want to apply an analytic description. If we assume 
the trapped UCN in mechanical equilibrium we can use formulas derived in the book [2j where 
the authors analyzed the effect of the earth's gravitational field on neutrons moving in a bottle. 
We adapt the notation to our case and replace the height parameter h by the radial coordinate r 
characterizing the multipolar magnetic field strength. At any point in the trap accessible for low 
field seeking neutrons with total energy eo the kinetic energy is given by E (r) = eo — \ V m (r)|. 
Note that we neglect the gravitational field, which is a good approximation for a horizontal 
source with less than 10 cm diameter. The angular averaged probability fx for UCN loss during 
a wall collision is a function of the kinetic energy E with which neutrons hit the wall. For a 
boundary between superfluid helium and the neutron optical potential V + iW with small losses, 
i.e. / = W/V <C 1, it is given by 



fx (E) = 2/ Re — arcsin J —-J — - 1 , V = V - U H e- (10) 




For convenience in later calculations we have included the projection on the real part of the 
expression to ensure that fx (E < 0) = without need to specify the range of E in advance as 
positive. The function ~p(E) rises monotonously with E from 71 = for E = to 71 = irf 
for E = V (for E > V we may set 7i = 1 since we are not interested here in calculating the 
dynamics of marginally trapped neutrons). For E = V'/2, 71 ~ 1.14/. Since we consider a long 
trap for which 2itRL » irR 2 , we can neglect losses due to W at the end disks, and the argument 
of 71 (E) becomes 

E = e - V mR . (11) 

The magnetic multipole suppresses wall losses because first, only a fraction of neutrons has 
sufficient energy to hit the lossy wall and second, because those with eo > V m R hit the wall due 
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0.159 


0.252 


0.465 


0.739 


0.824 


0.859 


0.89 


10 


0.248 


0.351 


0.555 


0.785 


0.855 


0.883 


0.909 


12 


0.325 


0.429 


0.619 


0.818 


0.876 


0.9 


0.922 


14 


0.39 


0.492 


0.667 


0.841 


0.892 


0.913 


0.932 



Table 1: Values for 7' as defined in eq. C[5]for various values of eo/Kn_R. 



to magnetic deceleration only with kinetic energy in the range < eo — V m R < Vt rap — V m R, and 
~p (eo — V m n) < "p (eo). A further suppression is due to a reduced rate of UCN wall collisions, 
which we calculate next. As a result of phase space transformation under the influence of a 
conservative potential, in the spatial region accessible for the low field seekers^, the spectral real 
space UCN density is given by 



n (e , t, r) = Re J £ ° |Fm(r)l n (e , t, 0) , (12) 
V e o 

where projection on the real part ensures n (eo, t, r) = for r > R* defined by eo = \V m (R*)\- 
We define an effective real space volume of the source for neutrons with total energy eo by 



7 (e ) = 2nL Re F J € ° |Kl (r)l rdr. (13) 
Jo V e ° 

The spectral UCN density averaged over the entire volume of the converter is then given by 

n(e ,t) = 7'(eo)n(e ,t,0), (14) 

and the reduced quantity 



V (eo) = ^ = 2 Re f Jl-^ri ~Vdr (15) 

was derived using eq. [H The values for 7' listed in table 1 show that the reduction of real 
space density with respect to a potential well of same depth becomes less significant for higher 
multipolarity. 

The spectral flux density of neutrons at any point in the vessel, per unit area and per energy 
interval about eo, is given by the gas kinetic relation 

J (e , t, r) = (e , t, r) v (e , r) . (16) 

2ttRLJ (eo, t, R) is the spectral rate of UCN collisions with the cylindrical wall of the helium 
container. The speed v(eo,R) of the neutrons as they hit the wall is related to the speed at 
r = through 

v (e , R) = Re J £ °~ VmR v (e , 0) . (17) 
V eo 



3 We do not consider the high field seekers further. They hit the wall with higher energies and will leave the 
trap much faster than the low field seekers. For a calculation of the polarization of the neutrons they should 
however be properly taken into account. 
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With eq. [12] we thus obtain 



J (e , t, R) = £ ° VmR J (e , t, 0) G (e - V mR ) , (18) 

where the step function (eo — V mR ) = 1 for eo > V mR , and otherwise. 

For the loss term in eq. [1] due to collisions with the cylindrical wall we can write 

" ( "'- /) -^Tl(e -V mR )J(e ,t,R). (19) 



T wall (eo) R 

Further evaluation using eq. [TSJ eq. [16] for r = 0, eq. [21 and inserting /Z from eq. [10] leads us to 



-l / \ t?(eo,0) F' / . /e - U mfl e -V mR f e - V mR \\ 

-wall M = -^T-/- Re V - \/ — — V - — — J J • (20) 

For the calculation of the corresponding expression for losses through slits it is reasonable 
to assume them to be situated at r = R, e.g. at the seam of the tube or at its connections to 
the circular windows for the cold beam. Assuming that any UCN hitting a slit will be lost and 
denoting the total surface of all slits by A their contribution to the loss term in eq. Q]is given by 

n(e ,t) 2 A 

7^)-R^RL J{e ^ Rh (21) 

neglecting the small surface of the disks at the ends. Hence, 

! Av(e ,0) e -V mR 
-slit (^) = 6 (6 - V mR ) , (22) 

and we see that the ordinary gas kinetic expression represented in the first fraction on the right 
hand side becomes reduced for V mR > for the same reason as the wall losses discussed before. 



4 UCN production 

In absence of the magnetic multipole field and for homogeneous irradiation of the converter with 
the cold neutron beam, the UCN production rate density is position independent, 

p =pl+pt= P (eo) de = KV*g. ( 23 ) 

Jo 

For an unpolarized beam, as always assumed hereafter, the production rates for the two com- 
ponents with spin parallel (t) and anti-parallel (J.) to the magnetic field are equal, i.e. = 
Pq = pq/2. The dependence follows for a homogeneous population of states within a sphere 
in momentum space. The factor K due to single phonon emission has been calculated on the 
basis of neutron scattering data and confirmed in several experiments [52], |66] [67], [68] , albeit 
with modest experimental accuracy limited by detection efficiency and other corrections. For 
unpolarized neutrons and UCN with maximum energy determined by Vtrap = V — Vn e = 233 
neV for Be or Ni with natural isotopic composition, it is given by 

K w 5s _1 cm" 3 $| 0-89mn [I0 9 cm _2 s~ 1 nm" 1 ] / (233neV) 3/2 , (24) 

where $o.89nm is the differential unpolarized neutron flux density at a neutron wave length of 
0.89 nm [69]. The flux unit is chosen numerically close to values available at existing facilities, 
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„ \ KniL 1 3 5 1 

" \ Vtrap 1 4 8 2 

4 0.229 0.37 0.456 0.551 

6 0.4 0.517 0.585 0.659 

8 0.512 0.609 0.665 0.725 

10 0.589 0.672 0.72 0.77 

12 0.645 0.718 0.759 0.802 

14 0.688 0.752 0.788 0.827 



Table 2: Values for k as denned in eq. [28] for various values of n and V^r/ Vtrap- 



e.g. the monochromatic beam H172a at the ILL [70]. An additional, usually smaller contribution 
to UCN production is due to multi-phonon processes [71] [72] . 

In the converter with the magnetic multipole switched on, the UCN production rate density 
becomes position dependent. For the low field seekers, 

p t (r) = * Re (y trap _ \ Vm ( r )|)3/2 (25) 

Re (...) accounts for the fact that in regions where \V m (r)\ > Vt rap no neutrons with total energy 
smaller than the trapping depth can be produced. After averaging over the volume of the 
converter vessel and using eq. [6] we obtain 

P t = *t£2 f Re f 1 " 'A rdr. (26) 

J V ^trap / 

Note that, with the definition of 7' (eo) given in eq. [T5l we can write the spectral mean UCN 
production rate density as 

pt( eo) = l*: 7 '( eo )^, (27) 

leading to eq. after integration over eo- 

The ratio of total production rates for low field seeking UCN with the magnetic multipole 
switched on and switched off, 

k = ^p, (28) 
Po 

is smaller than unity due to the phase space reduction by the magnetic multipole. The values 
quoted in table 2 demonstrate the positive effect of high multipole order n on k and hence on 
the saturated UCN density calculated in the next section. There are however practical limits. 
First, thermal insulation between the magnet and the much colder helium container necessitates 
an annular gap over which the field would drop too strongly if n is chosen too large. Second, 
the maximum field strength achievable with a given maximum current density in the current 
bars around the converter of given diameter decreases with n. For R = 5 cm, and taking into 
account the results given in the next section, n ~ 12 seems a reasonable choice, which also is 
still feasible using standard NbTi superconducting wire technology. 

5 Saturated UCN density 

The saturated spectral UCN density follows from eq. [2] with eq. [27] which already includes the 
necessary spatial average over the converter volume, 

nI (e ) = ^7 / (eo)v / ^r(6 ,T). (29) 
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Figure 2: Saturated density rioo of low field seeking UCN in a converter vessel with diameter 
10 cm, held at 0.5 K and surrounded by a 12-pole, respectively, 8-pole magnet. Calculations 
employ eq. [30] for Vtrap = V — Vn e = 233 neV (i.e. for a converter vessel made entirely of Be), and 
for various values of / = W/V . The solid lines show results for the best value of / previously 
achieved for Be, while the upmost curves show a hypothetical situation for illustration. Values 
are given for an unpolarized differential neutron flux density of ^o.SQnm = 10 9 cm _2 s _1 nm _1 at 
A = 0.89 nm, as available at the monochromatic beam position H172a at the ILL. A characteristic 
time constant To is calculated for neutrons with velocity v = ff m ax> using eq. [20] for V mR = 
and e = ^U tmp . 



Hence, using eq. 0and writing out all relevant arguments characterizing the multipolar magnetic 
field and the converter, the saturated total mean UCN density in the converter is given by 

[V mR ,n,VJ,V,T\ = -K / — — — v de . 30 

v ' 4 J t 1 (e , V mR ,n,V,f,T) 

The dependence on V is contained in the upper limit of integration, see eq. [7] From the various 
contributions to the rate constant r" 1 (see eq. [8]) we retain the terms due to wall collisions, 
upscattering and neutron beta decay, assuming that wall losses can entirely be described by eq. 
[20l and that there is no 3 He in the converter and there are no slits in the vessel. 

Figure 2 shows results for a vessel made entirely (i.e. V = V) of Be and exemplary for 
multipole magnets with multipole order n = 8 and n = 12, as a function of the magnetic 
potential at the cylindrical wall of the vessel, normalized to the trapping depth. The UCN 
density increases with n as expected due to the increase in trapping phase space. Beryllium 
has become a standard material for UCN trapping, with a best reported experimental value of 
/ = 3 x 10~ 5 in the low temperature limit [731 174|. despite much smaller theoretical values (the 
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Figure 3: Saturated density of low field seeking UCN in converter vessels with diameter 10 
cm, V = 252 neV, / = 3 x 10~ 5 , T = 0.5 K, and beam window and UCN valve (see fig. 1) 
coated with three different materials with values of V as indicated in the legend. The upper 
(lower) curve of each pair is for multipolarity n = 12 (n = 8), and all values are given for 
^o.89nm = 10 9 cm _2 s _1 nm _1 (unpolarized) . The kinks visible for the two upper pairs of curves 
appear at fields for which V m R = V — V. 

finding that these were never reached was termed "anomalous losses" and has triggered many 
experimental investigations and speculations). However, it might be more realistic to consider 
even larger values for /, assuming that efficient cleaning procedures cannot be applied in situ 
(e.g. baking is excluded in presence of indium seals). As obvious from the curves, the poorer the 
neutron optical UCN storage performance, the larger will be the improvement of UCN density 
due to the multipole magnet. For the situations reported in the introduction, with experimental 
loss rate ratios as high as 7" _1 /r^ 1 ~ 5.5 [34J, respectively 13 [39], a multipole magnet would 

indeed be very useful. Some of the curves exhibit a maximum of nto for values U m #/T4rap < 1- 
This can be understood as resulting from the competition of storage time constants r and the 
effective trap volume 7' entering eq. [301 For bad values of / the optimum obtains for V^ij/Vtrap 
close to 1, while for a trap with excellent storage properties the multipolar field reduces the 
UCN density even at low field values because the factor 7' < 1 then dominates over a marginal 
gain in r. For illustration of this behaviour we also added a curve for an unrealistic converter 
vessel with hypothetical / = 3 x 10 -6 . 

Next we consider an interesting further opportunity which takes advantage of the fact that 
the multipole magnet increases not only storage time constants but also the potential energy of 
the neutron at the cylindrical wall. Hence, the trapping depth of the converter vessel becomes 
larger if the disks providing axial confinement are made of a material with larger neutron optical 
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Figure 4: Saturated density of low field seeking UCN in converter vessels with diameter 10 
cm, V = 252 neV, T = 0.5 K, / = 3 x 10" 5 (solid curves) jar / = 10" 4 (dashed). Beam 
window and UCN valve (see fig. 1) are coated with c n BN (V = 315 neV). The vessels are 
surrounded by magnets with multipole order varied between n = 6 and 24. Values are given for 
^O.SQnm = 10 9 cm _2 s _1 nm _1 (unpolarized). 

potential V > V (remember eq.[7Jand see fig. 1). Since the surface of the disks is small, one may 
employ materials which would be unsuited for the entire vessel, for an unfavorably large value 
of / = W/V or because they might be difficult to produce or deposit in high quality on large 
surfaces. While diamondlike carbon has already been studied in some detail [75J E3 [77J [75], 
further candidate materials able to extend the spectrum for UCN trapping beyond the Be cutoff 
have been the scope of recent investigations [79] . Particularly promising is boron nitride in the 
cubic phase (cBN). Its neutron optical potential of 324 neV is even larger than that of diamond 
(304 neV) but due to the^ large absorption cross section of the isotope 10 B, cBN also has an 
excessively large value of/ = 1.5xl0~ 2 . Enrichment of the weakly absorbing 11 B however may 
reduce / down to 3.3 x 10 -5 , along with a further increase of V to a theoretic value of up to 350 
neV, depending also on the high-density cubic phase content of the specimen. Using experiments 
on transmission with time-of-flight analysis and cold neutron reflectometry, the authors of |79] 
have demonstrated that their 2 /im thick n BN deposit on a circular silicon waver possesses a 
value of 315 ± 10 neV. The deviation from the ideal value is mainly due to a cubic phase content 
of 90%, which was measured independently by IR spectroscopy. 

Figure 3 shows saturated UCN densities calculated for a Be trap with axial UCN confinement 
provided by Be, c n BN as reported, or cubic C3N4, the latter possessing an extraordinarily large 
theoretical value of V = 391 neV. Calculations were this time performed with a fixed value of 
/ = 3 x 10~ 5 . As a function of magnetic field at the cylindrical wall, the curves for V > V start 
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Figure 5: Saturated density of low field seeking UCN in converter vessels with diameter 10 cm, 

V = 252 neV, T = 0.5 K and surrounded by a 12-pole magnet. Solid curves are for / = 3 x 10 -5 , 
dashed ones for / = 10 -4 . V — Vn e is varied between 294 neV and 60 neV. Values are given for 
^o.89nm = 10 9 cm _2 s _1 nm _1 (unpolarized). 

with a slope larger than in the case V = V shown in fig. 2. This is due to the increase of Vt ra p 
and therewith of the integration range in eq. [30] for V m R < V — V. Kinks in the curves appear 
at V m R = V — V where the full trapping depth is reached. While for a high quality Be trap 
with V = V the gain in UCN density is not too impressive (lowest two curves), for a magnetic 
12-pole and V = 315 neV, a technically feasible field of Br = 2.5 T enhances the saturated low 
field seeker UCN density by a factor 2.2 from 820 cm -3 to 1820 cm -3 . The magnetic multipole 
is thus an attractive device even for a vessel with good storage properties. 

Figure 4 shows the dependence of saturated UCN density on the multipole order, for traps 
made of Be with / = 3 x 10~ 5 , respectively / = 10~ 4 , each with end windows made of c n BN 
with V = 315 neV. For the worse trap a field of Br = 2.5 T enhances the saturated low field 
seeker UCN density by a higher factor 3.6 from 390 cm~ 3 to 1400 cm~ 3 . Hence, while without 
multipole the worse / suppresses the UCN density by more than a factor two, with field it is 
only 22% lower than for the better trap. The multipole field thus stabilizes the output of the 
source; it becomes much less sensitive to the loss coefficient of the Be surface or if its quality 
degrades with time. 

Figure 5 shows, again for traps with the cylindrical section made of Be with / = 3 x 10~ 5 , 
respectively / = 10~ 4 , the dependence of saturated UCN density on the upper bound of the 
trapped UCN energy spectrum. This is easily calculated using eq. [30] by setting the potential 

V — Vue to appropriate values. We see that, with decreasing V, the magnetic field needed to 
optimize the UCN density decreases as well. The reason is that a converter with lower trapping 
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/ 3 x 1(T 5 1 x 1(T 4 2 x 1(T 4 4 x 1(T 4 

nlo (cm" 3 , for B R = 2.5 T) 1820 1400 1200 1040 

nlo (cm" 3 , without magnet) 820 390 230 130 

Table 3: Values for the saturated density of low field seeking UCN as defined in eq. [30] for a 
converter vessel with R = 5 cm, V = 252 neV, V = 315 neV, T = 0.5 K, surrounded by a 
12-pole magnet, and for several values of /. Values are given for an unpolarized differential flux 
density of <J?o.89nm = 10 9 cm _2 s _1 nm _1 . 



potential accumulates UCN with smaller energies for which wall collisional losses are smaller 
due to the energy dependence of Ji (E) defined in eq. [TO] and due to wall collisions occuring at 
a smaller rate. In addition, the effective volume 7' decreases with V m R faster than for neutrons 
with higher energies, as obvious from table 1. If one wants to extract (or the experiment 
connected to the source is able to use) only a spectrum with low-energy UCN of say, up to 60 
neV as provided for instance by a magnetic bottle with a trapping field of 1 T, the multipole 
magnet around the source will be useless for / = 3 x 10~ 5 . For the worse trap with / = 10 -4 it 
will however still lead to a gain factor attaining a maximum of 1.4 at Br ~ 0.75 T. 



6 Conclusions 

As our analysis shows, a multipole magnet can lead to a large gain in the saturated density of low 
field seeking UCN because the presence of the field reduces the number of neutrons hitting the 
material walls and reduces the energy and wall collision rate of those that do. A 12-pole magnet 
with field Br = 2.5 T on a radius of R = 5 cm seems technically feasible using existing NbTi 
superconducting wire technology, as investigated in an independent study using a finite element 
code. Based on results of experimental work done by other groups, a promissing candidate vessel 
consists of a Be trap closed off by disks coated with c n BN. Alternative materials are diamondlike 
carbon with V close to 300 neV depending on the abundance of sp3 chemical bonds, or enriched 
58 Ni with V = 346 neV and a theoretical / = 8.6 x 10~ 5 , which however is magnetic and UCN 
depolarization might be an issue that needs experimental study. We note that, in order to 
extract the full benefits, the incoming cold beam will need to be transported by a supermirror, 
with a top layer deposit of the UCN reflecting material with neutron optical potential V. An 
experimental investigation of this concept is currently underway at the ILL. 

Our benchmark converter is able to provide a saturated low field seeker UCN density almost 
as high as an unrealistic, perfect trap with / = and Be cutoff, for which one calculates 
nJo = 2060 cm -3 . In addition, the densities quoted in table 3 for various values of the parameter 
/ highlight a valuable practical advantage offered by a magnetic multipole reflector. Even for / 
by a factor 13 worse than the best value of 3 x 10 -5 previously achieved for Be, nto would still 
exceed 1000 cm~ 3 and only be a factor 1.8 smaller than in the best case. On the other hand, 
for a pure Be trap without magnet the same lack of performance in / would reduce the UCN 
density by a factor 6.3. Obviously, a converter vessel equipped with multipole magnet is then 
also much less sensitive to degradation of its surface quality with time. 

Up to now, we have considered only the low field seekers. Without the magnetic field the 
high field seekers are however equally well trapped, so that in this case (and still assuming 
the usual situation of an unpolarized cold neutron beam for UCN production) the total UCN 
density is a factor 2 higher. For experiments requiring polarized UCN such as a neutron lifetime 
experiment using magnetic trapping or the neutron EDM experiment, this factor is of no use 
and the benefit of the magnetic multipole mirror is obvious. 

As an important detail not affecting our conclusions we note that, in addition to the multipole 
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field it will be necessary to apply a weaker bias field in the order of some 10 mT along the 
converter axis to avoid depolarization in the region around r = 0, where the multipole field is 
zero. In addition, we can consider using axial magnetic pinch fields to increase Vtrap and thereby 
the density of the UCN, however for an extended energy spectrum. For extraction the field on 
one end needs to be ramped down, so that an iris type UCN valve might be more appropriate 
for this case. An extended UCN spectrum would be interesting if UCN of any velocity are 
beneficial for the experiment, such as UCN transmission experiments, or in combination with a 
phase space transformation by letting the UCN rise against the gravitational field. Note however 
that extraction of such a UCN spectrum will be a challenge. On the other hand, some studies 
might be performed in situ using static pinch fields, such as experiments on UCN upscattering 
in superfluid 4 He, for which any increase in UCN density will be very welcome. 
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